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Abstract 

A simple approach is proposed allowing actual calculations of the preequilibrium dynamics in 
ultrarelativistic heavy-ion collisions to be performed for a far-from-equilibrium initial state. The 
method is based on the phenomenological macroscopic equations that describe the relaxation dy- 
namics of the energy-momentum tensor and are motivated by Boltzmann kinetics in the relaxation- 
time approximation. It gives the possibility to match smoothly a nonthermal initial state to the 
hydrodynamics of the quark gluon plasma. The model contains two parameters, the duration of 
the prehydrodynamic stage and the initial value of the relaxation-time parameter, and allows one 
to assess the energy-momentum tensor at a supposed time of initialization of the hydrodynamics. 

PACS numbers: 25.75.-q, 24.10.Nz 
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I. INTRODUCTION 



A comprehensive analysis of the experimental data from the Relativistic Heavy-Ion Col- 
lider (RHIC) has shown that a quark-gluon plasma (QGP) is created in these collisions, 
and that the thermalized QGP is, perhaps, the most perfect liquid possible in nature. This 
conclusion is based, in particular, on the success of ideal hydrodynamics in describing of 
the basic features of heavy-ion collisions at RHIC energies (for a review see, e.g., Ref. 
Recently, essential progress has been made in the development and applications of viscous 
hydrodynamics for RHIC heavy- ion collisions (for reviews and recent results see, e.g., [l 3 1 
and references therein). Viscous hydrodynamics accounts for deviations from local equi- 
librium by means of dissipative transport coefficients, and, therefore, has a more extended 
region of applicability as compared to ideal hydrodynamics. Nevertheless, it should be em- 
phasized that the domain of validity of viscous hydrodynamics is still the hydrodynamic 
regime; thus hydrodynamics is valid when the relaxation time r rc i is much smaller than 
the inverse expansion rate, l/d^u^, that is, T Te \d^u^ <C 1, and one can expect that hydro- 
dynamics breaks down when r re i<9 M u M ~ 1/2. In practice, hydrodynamics breaks down in 
ultrarelativistic heavy-ion collisions at the very initial nonequilibrium stage, near the edge of 
the fireball, and at the later rarefied kinetic stage of matter evolution. Since hydrodynamics 
is based on the assumption that the system is near local thermal equilibrium ^-6|, it is 
natural that its region of applicability cannot be reliably determined from within itself, and 
can be properly estimated only with the help of an appropriate nonequilibrium theory. 

In addition to the breakdown times, one also needs to specify initial conditions, such as the 
energy density, fluid velocities, and viscous shear tensor to apply hydrodynamics. Evidently, 
the initial conditions for hydrodynamics are determined at the very initial nonequilibrium 
stage of the matter evolution in ultrarelativistic heavy-ion collisions. This stage is well 
understood now based on the color glass condensate (CGC) approach (for a review see, 
e.g., Ref. [zj]), which is adequate at RHIC and, probably, Large Hadron Collider (LHC) 
energies, although an explanation of the thermalization and nearly perfect fluidity provides 
a challenge for the theory of quark-gluon matter (for recent results, see Refs. [8|, |9(). It is 
worth noting, however, that the assumption of very early thermalization (say, r t h = 0.3 fm/c) 
is not necessary for data description in hydrodynamics, since the transverse collective flows 
and their azimuthal anisotropy for noncentral collisions can appear already at the beginning 



2 



of the hydro dynamic expansion as a result of the development of the transverse velocities at 



the prethermal (glasma 10] or partonic or string stage 12] . That is why an analysis 
of matter evolution at the prethermal stage, which determines the further hydrodynamic 
expansion, is so important.. 

The nonthermal initial conditions are related typically to the Bjorken proper time Tq 
when the system can be characterized by the phase-space density of individual partons. 
In the CGC approach this "formation time" was estimated to be To ~ 0.1 — 0.3 fm/c for 
central collisions in midrapidity at LHC and RHIC energies [l3, [ijj]. Taking into account 
the theoretical estimate of the thermalization time scale r t h = 1 — 1.5 fm/c for the LHC 
and RHIC heavy- ion collisions [9j, one can hardly expect hydrodynamics to be applicable 
at r = t < 1 fm/c, because then r re id^u^ ~ T th /r > 1. 

Thus, to obtain the initial conditions for hydrodynamics at r = r t h > 1 fm/c, one needs 
to match the very early initial stage of the nuclei collisions, when the hydrodynamic approx- 
imation is outside the regime of its validity, with an almost local equilibrium state of the 
QGP. It is noteworthy that phase-space distributions derived from the CGC approach are 
highly anisotropic, and so the distribution functions are quite nonequilibrium. The longitu- 
dinal momentum distribution is much more narrow than transverse one, and asymptotically, 



at late times, it is a delta function 6(p z ) at z = [15j. The latter approximation is often 
used for a description of the initial stage in A + A collisions (see, e.g., 16J). 1 

In this paper, we formulate phenomenological macroscopic equations accounting for the 
energy-momentum and quantum number conservation laws which allow one to connect the 
arbitrary initial, probably highly nonequilibrium state of the system with the (partially) 
equilibrated thermal state, when one can start to use (viscous) hydrodynamics to describe the 
further matter evolution. 2 The model exploits a minimal set of phenomenological parameters 
and can be utilized to assess the initial conditions for hydrodynamic expansion, which can 



1 One should understand, however, that such a utilization of the asymptotic approximation is not quite 
correct from the point of view of quantum theory. Indeed, the Wigner function fyj(x,p), which is 
the quantum mechanical analog of the classical phase-space density f(x,p), satisfies the restriction 
J f^(x,p)d 3 pd 3 x < (2ith)~ 3 (see, e.g., nere the normalization condition J fy/(x,p)d 3 pd 3 x = 1 
is supposed. Then, in order to escape a contradiction with quantum mechanics, smooth boost-invariant 
prescriptions for the longitudinal part of the distribution f{x,p) from Refs. [18] or 19] should be used 

for a description of the initial stage in A 4- A collisions. 

2 A similar effort was undertaken in Rcf. 



20(, where the preequilibrium dynamics of a (0 + l)-dimcnsional 



QGP was matched to second-order viscous hydrodynamics. 
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be used, then, for a hydro dynamic modeling of ultrarelativistic heavy-ion collisions. 



II. ENERGY AND MOMENTUM RELAXATION DYNAMICS FOR THE EARLY 
STAGE OF ULTRARELATIVISTIC HEAVY-ION COLLISIONS 

It is well known that statistical derivation of the macroscopic hydrodynamic equations 
for gases is based on an approximate solution of the Boltzmann equations near the local 



equilibrium distribution (see, e.g., Refs. The form of hydro* lynamio equal ions derived 

can then be spread out phenomenologically from a gas to dense fluids if they are near 
local thermal equilibrium. More precisely, hydrodynamic equations for dense systems like 
liquids near local thermal equilibrium can be derived by means of the (quasi) equilibrium 
statistical operator method based on the conditions of a maximum of local entropy and 



under appropriate auxiliary conditions js]. The proximity to local thermal equilibrium 
allows one to utilize for a description of these inhomogeneous evolutional systems the Gibbs 
thermodynamic relations in the same form as for equilibrium systems [5j, and to get the 
equations of viscous hydrodynamics in the closed form. However, if the system is far from 
local equilibrium, it can be described macroscopically only through the unclosed hierarchical 
structure of the hydrodynamic balance equations expressing the conservation laws in the 
system 5|. All the momenta of the physical values are present there, and they cannot be 
calculated without knowledge of the evolution of the nonequilibrium distribution function 
4j . Although the second-order viscous hydrodynamic equations increase precision compared 
with the first-order equations and, moreover, have attractive features from the point of 
view of mathematical formulation, the domain of applicability of these equations is still the 
hydrodynamic domain: the microscopic scales like the mean free path are much smaller than 
the macroscopic dimensions (the system's homogeneity lengths). Because of this, one can 
conjecture that dissipative fluid dynamics is inapplicable in its standard form to the very 
initial stage, r < 1 fm, of ultrarelativistic heavy-ion collisions (see also Ref. when 
the rate of expansion is too high and the interactions are not strong enough to validate 
the application of the hydrodynamic approach. Thereby, the very initial nonthermal state 
(glasma?) of matter evolution in A + A collisions cannot be described hydrodynamically, 
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unlike the subsequent thermal QGP stage, where hydrodynamics seems to work well [l|. 3 

To perform the heuristic arguments to derive phenomenological macroscopic equations 
for a description of transient prehydrodynamic behavior, let us first assume that the initial 
phase-space density is associated with the partons, while the resulting expressions can be 
used as a phenomenological approach not only for partonic matter but also for other forms 
of prethermal matter like glasma and strings. Usually, if the exact dynamics is cumbersome 
or unknown, establishment of the local equilibrium state in the system is modeled by the 
kinetic equation for the phase-space distribution function f(x,p) in the relaxation time 
approximation, with f leq (x,p) as the target function to which the phase-space distribution 
tends. The corresponding equation has the form 

p>*df(x,p) _ M f(x,p) - f cq (x,p) 

where r* el is the relaxation-time parameter in the local rest frame of the energy flow (in 
general it is some function of (x,p)), and u^(x) is the four- vector energy flow field, 

T^(x)= [ d 3 p^f(x,p). (3) 
J Po 

To determine the parameters (temperature, etc.) that define the local equilibrium state, 
f lcq (x,p), it is necessary to implement the energy-momentum conservation equations, 

d,T^{x) = 0, (4) 

and, if necessary, the conservation equations for (net) quantum numbers qi, 

d^(x) = 0, (5) 

with current qf(x), 

^k»[f qi (x,k)-f_ qi (x,k)}. (6) 



3 Note that the results of viscous hydrodynamic calculations depend on exact values of the transport coef- 
ficients that still cannot be derived unambiguously from the complicated theory of microscopic dynamics. 
At the later, low-density gaslike stage, the viscous hydrodynamics can be coupled to the Boltzmann 



kinetics of the hadron resonance gas 22j | . 
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Then Eqs. (J4]) and (jSJ) together with the formal solution of Eq. ([I]) , 

f(t, r, p) = /(t , r-^-(t- t ),p)P{t , t, r, p) + 
Po 



E 

/ /tr"-^-^?)^^^)^, (7) 
J po at' 

to 

where the probability for particle with momentum p to propagate freely from point (t', r) 
to point (t,r + jp(t — t')) is 

{V p u uHs,T- *-(t- s)) ) 
J Por* el (s,r-^(t-s),p) J 

lead to very complicated equations for the thermal and hydrodynamic fields (tempera- 
ture, collective velocities, etc.) that define the target local equilibrium state described by 
f leq (x,p). The problem becomes even more severe if one takes into account that r* el is a 
function of /' eq . Also note that the target state is reached in a finite time interval at t = t t h, 4 
only if the relaxation-time parameter in Eq. ([I]) vanishes at t — > t t h- T* el (t — > t t h, r,p) — > 0. 

All this makes difficult a utilization of Eq. ([1]) for a matching of a far-from-equilibrium 
initial state to perfect or viscous hydrodynamics with presupposed transport coefficients 5 in 
relativistic heavy-ion collisions. The situation reflects the fact that the solutions of the Boltz- 
mann equations as a rule do not correspond to any hydrodynamics, or, in some conditions 
(in the vicinity of a local equilibrium state), correspond to solutions of viscous hydrody- 
namic equations with a very specific form of the viscosity coefficients defined by the cross 
sections or by the relaxation-time parameters in the relaxation-time approximation of the 
collision terms. To utilize this in practice a quite complicated Chapman-Enskog method 



231 ]) is applied. This method is inapplicable, however, if the system is initially 



(see, e.g., 
nonthermal. 

Here we propose a simple phenomenological method allowing one to describe the relax- 
ation of an initially far-from-equilibrium system toward the neighborhood of a local equilib- 
rium state. The latter is described by the target energy-momentum tensor and is associated, 



4 Of course, the same statement is easy to reformulate for the case when one operates with an arbitrary 
three-dimensional spacelike " thermalization" hypersurface where t — tth(r) instead of a constant time t t h- 

5 Then the target function in Eq. (jTJ) could have a nonequilibrium form with parameters that are linked to 
transport coefficients. 
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in general, with viscous hydrodynamics. For simplicity, we will consider in detail a model 
with the target energy-momentum tensor associated with a perfect fluid; a generalization for 
the viscous target energy-momentum tensor is straightforward. To arrive at such a model, 
it is useful to start again with the relaxation-time approximation ([I]) which has the formal 
solution (JZJ) , (E]) - Then, assuming that the maximum of the integrand in ([7]) occurs at the 
upper limit of the integral and, moreover, that / lcq (x,p) is a relatively smooth function, one 
can factor out from the integral in ([7]) the function f leq (t',r — ^(t — t'),p) at point t' = t 
and arrive at the simple form, 

f(x,p) = f(t ,r - ^-(t - t Q ),p)P(t ,t,r,p) + / leq (t,r,p)(l - P(t ,t,r,p)), (9) 
Po 

which demonstrates an explicit transition to a local equilibrium state at t = t t h if 
P(t , t th , r,p) = and is, of course, just an approximation of (j7j). Let us write this expression 
in some more general relativistic invariant form, aiming to use the curvature coordinates in 
our further analysis, 

f{x,p) = fi^{x,p)V{x,p) + f lc<1 (x,p)(l -V(x,p)), (10) 

where fi Te e(x,p) is the distribution function of freely streaming partons, 

P^ M /free(z,p)=0, (11) 

which coincides with f(x,p) at some three-dimensional spacelike hypersurface where initial 
conditions for f(x,p) are formulated, for example, ff Tee (x,p) = /(to, r ~ &(t — to),p) for the 
initial distribution function f(t ,r,p). Here V(x,p) has the same probability interpretation 
as P: it is the probability for the particle emitted on the initial hypersurface (in simplest 
case - at time to) with momentum p to reach point x: 

P{x,p) = exp{-J ^(x-x')M|J, (12) 
p»d„G p (x - x') = 5 (4 \x-x'), (13) 

V = 1 initially at t = to (or at some three-dimensional hypersurface to( r ) where initial 
conditions are specified), V — at t — t t h (or at some hypersurface tth( r )) when r* cl — > 0. 
One can easily see that (|T0|) satisfies the equation 

p»df(x } p) _ f(x,p) - f eq (x,p) p»df Xc \x,p) 

dx» ' P M r* el (x,p) dx» 1 1 ,P)) ' 1 J 



Note here that, while Eq. (TTOj) is relatively simple, the equations for f leq (x,p) that follow 
from conservation laws (j3J) and (j3J) are, in general, very complicated (even if one does not 
relate r* el with f leq (x,p)) because of the momentum dependence of V(x,p). 

Our idea is to utilize the approximate formal solution (fTTjj) . which preserves the most 
important properties of the true dynamics at the prethermal stage of evolution, to describe 
the relaxation dynamics of the energy- momentum tensor. To do it, we first make the ap- 
proximations 

r T * cl (x,p) « r r * cl (x), V(x,p) « V(x), (15) 

which significantly simplifies modeling of the relaxation dynamics at the early prethermal 
stage of evolution in ultrarelativistic heavy-ion collisions, 6 because then, in particular, the 
energy-momentum tensor (j3j) takes the simple form 

T^(x) = TgJx)V(x) + T&(x)(l - V(x)), (16) 

T£Ltx)= I d 3 p^f ircc (^ P ), (17) 
J Po 

T;; d (x) = / d 3 p^f^(x, P ). (is) 

J Po 

Even though (fl6|) is a rather rough approximation for the true energy-momentum tensor, it 
preserves the desired properties of the true expression, namely, it demonstrates a transition 
from a far-from-equilibrium state to a local equilibrium state and, which is very important, it 
allows one to account for the energymomentum conservation laws and, if necessary, quantum 
number conservation laws in a simple form. That is, accounting for the energy-momentum 
conservation laws <9 M T Miy (x) = 0, with T^ v {x) defined in (Tl6l) . and taking into account that 
d ll T^ e (x) = 0, we arrive to the following equations: 7 

0„[(1 - V{x))T£; d (x)\ = -Tr cc (x)d,V(x). (19) 

6 The example of explicit representation of V(x) will be given in the next section in hyperbolical Bjorken 
coordinates for constant longitudinal proper times corresponding to initial and final (" thermalization" ) 
hypersurfaces. 

7 Note here, to avoid misunderstanding, that while the expression (TT6|) looks like a smooth interpolation 
between the free-streaming and hydrodynamic regimes, it does not mean that initially there are no inter- 
actions and then interactions gradually switch on. One can see that this is not the case from the fact that 
the expression (flT)]) is based on the approximate formal solution (TlT)]) of kinetic equation (fl]) . 
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The initial conditions for T^ d follow from the right-hand side of Eq. ( fl4|) . That is, taking 
into account approximation ffl5|) and utilizing the fact that initially V = 1, we see that the 
initial conditions are given by the equation 

T^(x in )u^x in ) = T^ d (x in )u^x in ), (20) 

where x- in = (to(r),r). Note that the same equations defining the initial conditions for T^ d 
follow from Eq. ([1]). The solution of Eq. (120]) is straightforward and follows from the 
definition fl2]): 

ChydOin) = e(x in ), K yd (x in ) = ^(Xin) (21) 

As for the hydrodynamic pressure and equation of state (EOS), phyd = Ph y d( e h y d), they 
can be chosen, for example, in agreement with the lattice QCD EOS for the RHIC and 
LHC energies, and the same for the viscous coefficients. Then one can use Eq. (Tl9~]) with 
the initial conditions ( l2Tj) as phenomenological macroscopic equations for a description of 
transient pre-hydro dynamic behavior in ultrarelativistic heavy-ion collisions. 



III. A SIMPLE MODEL TO ASSESS THE INITIAL CONDITIONS FOR HYDRO- 
DYNAMICS 

Let us begin with the example, which we choose merely because of its simplicity, where 
the initial and final conditions are specified at constant-time hypersurfaces t = t Q and t — t t h, 
respectively. Then, in obvious notations, 

„ 9 d d /nnS 

p ^ = p »a +p a? (22) 

and we get from ( fl3i) and ( 1221) that 

G p {x-x') =j9 o 1 0(t-t , )©(i'-io)5 (3) (r(t,t / ) -r'), (23) 
where 0(t' — to) indicates that the evolution time starts at to, and 

T(t,t') = T-(p/ P0 )(t-t') (24) 

satisfies the equation 

p"d„r(t,f) = 0. (25) 
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Substituting these results into (fT2l) and integrating over r' we get an expression that explicitly 
coincides with ([8]): 

nt,T,p) = ev{-'[ * (26) 



to 



where 



irel(s ,r(M). P ) = p f' ( ;- r( ;; 3) :f . (27) 

^w M (s,r(t,s)) 
Referring to Eqs. (fTTj). (j22|) and (123]). we see that 

fbee{x,p) = f(t ,r(t,t ),p). (28) 

Equation (1281) allows to build the evolution of T^ e and after specification of r* cl to fix the 
equation of relaxation ( fl9|l . In the next example, which is very important in practice, we 
analyze this procedure in detail. 

Now let us consider the hypersurfaces r = r and r = r t h as the initial and final ("ther- 
malization" ) hypersurfaces, respectively; here r = y/t 2 — z 2 is the longitudinal proper time. 
In order to describe the boost-invariant dynamics, it is convenient to switch to the coor- 
dinates (r, rT,r]), where is thetransverse radius vector and r\ is the space-time rapidity, 
i] = tanh." 1 z/t; then t = rcoshr/ and z = rsinhr/. The particle four-momentum can be 
expressed through the momentum rapidity y = tanh^ 1 pi/p where p^ is the longitudinal 



momentum, the transverse momentum p^, and the transverse mass = ^Jm 2 + p\, then 
_ (j JlT cos hy ; p T) rrtT sinhy). Accounting for the boost-invariant dynamics, we choose 

T rei(x,p) = r* el {r, 6, r T , p T ), (29) 

where 9 = y — rj. Transforming p^d^ to these coordinates, we find 

P^tt- = mr cosh 0— H sinhfl— + p T - — . (30) 

ox 11 or t or] or T 

The auxiliary function G p then reads 



&(r-r')&(r'-r )5(r](r,T')-r]')6^(r T (r,r') 

cosh(?/ — r/J 

where 0(r' — r ) indicates that the time evolution starts at nonzero proper time r . Here 
r](r,r') and r T (r, r') satisfy the equations 

2AVKr,r') = 0, (32) 
^«V- T (r,r') = 0, (33) 
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and can be written as 

sinh 9{r, r') = — t sinh 9, (34) 

T 



r T (r,r') = r T - — (r coshfl - Vt' 2 + r 2 sinh 2 0), (35) 
my 

where 9(r, r') = y — r](r, r'). Substituting fl3Tl) into f|T2|) . accounting for (1291) . and integrating 
over v( and r^, we get 

V{t, 9, r T , pr) = exp - / - ^— -dr' , (36) 

^ J T lA {T , ,9(T,T , ),T T (T,T i ),p T ) J 

where 

r rcl (r',0(r,r'),r T (r,r'), P T) = - pTl^f^ ^ T r * el (-r', g(r, r') ; r T (r, T'), Pt)- (37) 

P^(t',9(t, r'),r r (r, r')) 

Then f iTee (x,p) is given by 

fhee(x,p) = /(ro,6»(r,r ),r T (r,ro),pT). (38) 

In order to pass now to a description of the relaxation dynamics at the prethermal stage, 
let us take into account that at this stage the longitudinal flow is much stronger than the 
transverse one and so, according to Eq. (1371) . r re \ m r* el . Therefore, to get a convenient 
parametrization of the prethermal evolution of the energy-momentum tensor of the system, 
one might assume that the relaxation time r* el in the rest frames of the fluid elements 
depends mainly only on the proper time r: r* el = T* el (r). Then V also depends on r only. 



It corresponds to the Bjorken picture 24jj where the thermalization processes are supposed 
to be synchronous in proper time of the fluid elements, so that the complete thermalization 
and the beginning of the hydrodynamic expansion happen at some common proper time Tth- 
Within such an approximation we get 

T^(x) = TZ(x)V(t) + T^ d (x)(l - V(t)), (39) 

where 

TO*) = J d 2 p T d9p^f(r , 9(r, t ), r T (r, r ), p T ), (40) 

and Tj^j(x) is associated with ideal or viscous hydrodynamics; for example, for the former 
it can be parametrized as follows: 

T hyd( x ) = (ehyd(x) + p hyd (x)K yd (xK yd (x) - p hyd (x)g^, (41) 
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where eh y d(^) is the energy density in the comoving system and Ph y d(x) is the pressure. The 
corresponding evolutional equations are 8 

- V{r))T^{x)\ = -T { Z(x)d,V(r). (42) 

In order to link Eqs. fl39|) and f|42l) with the hydrodynamics at r = T t h, one needs 

V(r ) = l, V(r th )=0, d fl V(Tt h ) = 0. (43) 

The initial conditions for T f ^ e coincide with initial conditions for T^, namely 
^free( r 05 t t, v) = T^" (tq, tt , rf) . The initial conditions for T^ d (see Eq. (120]) ) are defined 
by Eq. fpZTl) . The hydrodynamic pressure and equation of state, ph y d = Ph y d(eh y d), can be 
chosen, for example, in agreement with the lattice QCD EOS. 

In order to apply the equation of relaxation dynamics (T4"2"l) to calculate the space-time 
evolution of the energy-momentum tensor toward the hydrodynamic one, the function V(t) 
has to be specified. Without a knowledge of the specific thermalization dynamics, we propose 
to make the approximation for this as simple as possible with a minimal set of parameters. 
With that end in view, we do not discuss here the most general case, but propose the 
following simple ansatz for Vir): 



V(r) = exp{- I ^ds}, Hi) 



TO 



where 



Trei(s) = r re i(r )— — -. (45) 
7th - T o 



Performing the integral in ff44l) . we find that 

T th~ T 



Tth - TO 



V{t) = ^— . (46) 



Note that Eq. (|42p can be considered in computational dynamics as the hydrodynamic equation for the 
energy-momentum tensor = (1 — 'P(t))TjJ 1 d of an ideal fluid with an explicit "source" term on the 
right-hand side and with rescaled energy density £hyd = (1 — 'P(T"))ehyd, and pressure phyd = (1 — ^(^Phyd- 
Also, since in our approximation the relaxation time and, hence, the probability V of freely propagation 
are supposed to be known functions (up to the two fitting parameters), as well as the tensor of free 
propagating system T^ e (x), then Eqs. (j42|) are of the same type as the equations d^T^^ix) — 0. So, if 
the latter are equations of the hyperbolic type, as in the case of the ideal fluid or second-order viscous 
hydrodynamics (in the form of Israel and Stewart), then causality is preserved in our approach. 
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The conditions ( 143|) require > 1 (then also V/r* el — > at t — > r th ) which is a constraint 

on the model parameters: the lifetime of the prehydrodynamical period (r t h — t ) and the 
initial value of the relaxation-time parameter r re i(ro). The energy-momentum relaxation 
equation ( H2|) together with the proposed parametrization f j4~6l) for the probability function 
"P(r) can be used to assess T MV (r t h, iy, r/) = T^ d (rth, ?t, v) an d so find the initial conditions 
(in particular, the energy density and transverse flow) for hydrodynamic simulations of 
ultrarelativistic heavy-ion collisions. 

The heuristic method, which allows us to derive these macroscopic equations, is based 
on the Boltzmann equation in the relaxation-time approximation, but, to our mind, its ap- 
plicability area is not restricted to a rarefied system just as the form of the hydrodynamic 
equations based on approximate solutions of the Boltzmann equations near the local equi- 
librium distribution can be applied for a description of dense fluids if they are near a local 
thermal equilibrium. It is worth noting that Eqs. fj4~2]) . ( 144")) . and (1451) do not contain de- 
tails of the system's properties, except for the initial energy-momentum tensor and its free 
evolution. The method makes it possible to match the initial nonequilibrium state of the 
system and the locally equilibrated one for various kinds of system, such as partons, strings, 
and a glasma field. It also allows one to match a far-from-equilibrium initial state not only 
to the locally equilibrated one but also to the state associated with the initial conditions 
for viscous hydrodynamic evolution of the QGP. For the latter, Eqs. fj4~2]) and f|46|) preserve 
their form, but the energy- momentum tensor T^ d (x) differs from the simple form of perfect 
hydrodynamics ( l4Tj) and corresponds to the tensor of viscous hydrodynamics. It should be 
noted that in this case Eq. (1431) loses its interpretation as the relaxation time to the local 
equilibrium state (it tends to zero at r — > r th ) and so the probability (1461) plays now just the 
role of the interpolating function. Nevertheless, it allows one to match smoothly the very 
initial state of the matter in A + A collisions with the initial conditions for viscous hydro- 
dynamics when (partial) thermalization is already established; in other words, to find the 
initial hydrodynamic parameters - energy density, hydrodynamic velocity field, and so on - 
for the dissipative evolution of the QGP in agreement with the model relaxation dynamics 
and the conservation laws utilized at the prethermal stage. 
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IV. SUMMARY 



In this paper we presented a simple model of the early stage in ultrarelativistic heavy- 
ion collisions. The model describes smooth conversion of the preequilibrium dynamics into 
hydrodynamics, either perfect or viscous. Our phenomenological approach is motivated by 
the Boltzmann equation in the relaxation-time approximation, accounts for the energy and 
momentum conservation laws, and contains two parameters: the lifetime of the prehydro- 
dynamic stage and the initial value of the relaxation-time parameter. The preequilibrium 
evolution is modeled by the continuous evolution of the energymomentum tensor from the 
initial far-from-equilibrium state to the perfect or viscous fluid form. Themodel allows the 
flows and energy densities to be assessed at a supposed time of initialization of hydrody- 
namics, which then can be used as the initial condition for hydrodynamic simulations of the 
further evolution of matter (QGP) in ultrarelativistic heavy-ion collisions. 
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